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Abstract 

Using a reduced model focusing on the in-plane dependence of plane Couette flow, 
it is shown that the turbulent — > laminar relaxation process can be understood as 
a nucleation problem similar to that occurring at a thermodynamic first-order phase 
transition. The approach, apt to deal with the large extension of the system considered, 
challenges the current interpretation in terms of chaotic transients typical of temporal 
chaos. The study of the distribution of the sizes of laminar domains embedded in 
turbulent flow proves that an abrupt transition from sustained spatiotemporal chaos to 
laminar flow can take place at some given value of the Reynolds number i?i ow , whether 
or not the local chaos lifetime, as envisioned within low-dimensional dynamical systems 
theory, diverges at finite R beyond R\ ow . 

PACS: 47.20.Ft, 47.27.Cn, 05.45.Jn 

The transition to turbulence in flows lacking linear instability modes, such as Poiseuille 
pipe flow driven by a pressure gradient along a circular tube (Ppf) and plane Couette flow 
driven by two plates moving parallel to each other in opposite directions (pCf), is particularly 
delicate to understand owing to its abrupt character, without the usual cascade seen in the 
globally super- critical for e.g. convection. A recent general presentation of the issues 

is given in pQ. These globally sub-critical flows become turbulent through the nucleation and 
growth or decay of turbulent domains called puffs (Ppf) or spots (pCf), see, e.g., [2] and [3] 
for Ppf and pCf, respectively. Most of the work on the transition problem has dealt with 
special nonlinear solutions (exact coherent structures [U b]) to the Navier-Stokes equations 
and their dynamics in the phase space spanned by them. Such solutions, obtained within 
the so-called Minimal Flow Unit (MFU) assumption [5], have been found at moderate values 
of the Reynolds number R in Ppf, [6], pCf jHEHS], and also in plane Poiseuille flow [8|fT0]. 
(In pCf, R is defined as Uh/v where U is the speed of the plates driving the flow, 2h the 
gap petween the plates, and v the kinematic viscosity of the fluid.) These solutions are all 
unstable and, together with their stable and unstable manifolds, they form the skeleton of 
the turbulent flow at a local scale. The reason why turbulence can be sustained only at 
much higher values of R (about a factor of two to three higher) is however not clear [2]. 
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In practice, the existence of these nontrivial solutions has mainly served to explain the 
exponential behavior of the distribution of lifetimes of transient turbulence in the low-i? part 
of the transitional regime in terms of unstable periodic orbits, homoclinic tangles and chaotic 
repellers [11]. Unfortunately, this does not bring a definitive answer to the controversial 
issue of the existence of a threshold i?i ow above which turbulence could be sustained, see 
[2l[T2l[T3] for Ppf and [T2"j ll~i] for pCf. The reason is that the approach is zero-dimensional 
in essence (dynamics is condensed onto a small set of nonlinear modes with well defined 
spatial structure), which would be appropriate if confinement at the size of the MFU were 
really effective [151 §2.3]. In fact, the transitional flows considered are open and, ideally, 
unbounded in one spatial direction (Ppf) or two (pCf), so that a genuinely spatiotemporal 
dynamics is expected. 

A different angle of attack was proposed long ago by Pomeau [T6] (also [T71 Chap. 5]) who 
indeed stressed the importance of the spatial extension of systems experiencing subcritical 
bifurcations. The coexistence of states in phase space, typical of subcriticality, translates 
into coexistence of states in physical space, with walls separating homogeneous domains of 
each kind. He also pointed out that, when one of the competing states was chaotic and 
transient, one could expect a stochastic wall propagation similar to directed percolation |18j . 
another contamination process, with the obvious identification: locally turbulent — > active 
(alive), locally laminar — > absorbing (dead). According to this scenario, called spatiotemporal 
intermittency (STI) by Kaneko [19], transient chaotic local states of a distributed system 
may evolve into a sustained turbulent global regime due to spatial coupling. STI has been 
much studied within the framework of critical phenomena in statistical physics, especially 
in view of its universality [20J. Coupled map lattices in different spatial dimensions were 
used to test STI's properties [2"Tll2"2"] . In two dimensions, such an abstract minimal model 
was found to display a discontinuous first-order like transition as the control parameter 
was varied [23J, which was subsequently used for comparison with transitional pCf [HI b]. 
Unfortunately, it was too far removed from hydrodynamics to be really relevant, which 
prompted the derivation of a more realistic model directly from the Navier-Stokes equations, 
to be used here for studying statistical properties of the turbulent — > laminar transition in 
domains of very large lateral extension. 

A system of partial differential equations was derived by means of standard Galerkin 
expansion/projection of the Navier-Stokes equations, using a polynomial basis appropriate 
to no-slip boundary conditions; see [21] for a full description and details on its numerical 
implementation. This system extends Waleffe's model [4] by unfreezing the in-plane spatial 
dependence of the velocity field while keeping a few wall-normal modes. When restricted 
to lowest significant order, as is done for results discussed here, the model involves three 
fields. The first one describes streaks that are velocity perturbations pushing and pulling 
vortical structures. The two others account for the rest of the flow, especially the streamwise 
vortices which, together with the streaks, are believed to form the main ingredients of the 
turbulence self-sustaining process [US]. The model displays the most essential properties of 
the Navier-Stokes equations including lift-up, energy-preserving advection, and dissipation; 
it is also readily shown that the base flow remains linearly stable for all Reynolds numbers. 

When the streamwise length L x and spanwise width L z of the cell are large, simulations 
closely mimic the behavior of pCf in extended geometry. A direct transition toward a 
homogenous turbulent state is observed at sufficiently large R, characterized by a turbulent 
energy density independent of L x , L z > 16 (lengths given in units of the half gap h) and 
by a uniform streamwise velocity correction giving its typical S shape to the turbulent 
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mean flow [2U a]. This essential correction is also present inside and around developing 
spots J25J b] as observed in the experiments, see also |25J b]. However, due to the model's 
limited cross-stream resolution which underestimates energy transfer to, and dissipation in, 
small-scale wall-normal fluctuations, the transitional regime observed in the model takes 
place at R ~ 170, compared to R ~ 325 in laboratory experiments. 

In the following, the discussion is limited to statistical properties of the model at the 
turbulent — > laminar transition as obtained from numerical simulations, using a conventional 
Fourier, fully de-aliased, pseudo-spectral code second order in time, appropriate for periodic 
in-plane boundary conditions. Various spatiotemporal resolutions were used with similar 
results provided that numerical stability was guaranteed. (A subcritical bifurcation toward 
a high-energy numerical mode was observed when the spatial step was too large.) Domains 
of various aspect ratios ranging from moderate to large, and then very large, have been 
considered. All of our domains are much larger than the size of most active small scale 
structures in the system, which is ~ 6 x 3, the MFU. The size of the experiments [T3J,[2n] to 
which we compare our results ranged from 380 x 70 to 770 x 340. 

Moderate aspect ratios, 32 x 32, were first considered in |24, a], showing turbulent tran- 
sients with exponentially distributed lifetime probabilities. The characteristic times of these 
distributions were observed to increase exponentially as in some Ppf experiments [12J but 
seemingly faster for R approaching 175, suggesting a possible divergence around this value, 
i.e., the existence of a real threshold i?i ow ~ 175. There were indications that this apparent 
cross-over could be due to spatiotemporal effects. 

A larger domain, 256 x 128 was next considered [27]. It was first observed that the 
turbulent state could be maintained without difficulty down to R ~ 170 over time intervals 
at least of the order of 3 x 10 4 advective time units, and that the average total perturbation 
energy per unit surface E t fluctuated wildly [27, Fig. 7, left], with highs corresponding to 
a homogeneous distribution of small laminar patches immersed in turbulent flows and lows 
to a much more inhomogeneous situation with a few wide laminar domains surrounded by 
the same fine-grained mixture as for highs. Fluctuations of E t could then be attributed to 
the growth and decay of wide laminar domains. In [2U a] it was shown that the turbulent 
state maintained itself in the 32 x 32 system as long as E t (t) stayed above ~ 0.02 (in 
units of U 2 per unit surface), while it inevitably decayed below. Histograms of E t (t) were 
constructed for 32 transients observed for R = 170, Fig. QJa), and compared to those of 32 
time series of E t obtained by averaging over 32 sub-domains of size 32 x 32 covering the 
(256 = 8 x 32) x (128 = 4 x 32) system at the same R, Fig. Q](b). The histograms of the 
transient time series are less well defined than those of the sub-domains time series because 
the transients were always much shorter than the experiment in the 256 x 128 system. All 
show similar humped shapes and identical behavior beyond their maxima but at, and below 
E t « 0.02, an accumulation of values can be observed for the sub-domain time series whereas 
the histograms drop to zero for the transients. Since the laminar patches observed in the 
256 x 128 system were much wider than 32 x 32, the interpretation is obvious: a wide domain 
of size larger than 32 x 32 could become temporarily laminar (E t < 0.02) but was bound 
to return to the turbulent regime by contamination from its neighborhood while the small 
system could not recover (the transient ended), hence the link to directed percolation [16J. 

The first-order features of the transition have been examined by performing simulations 
in an even wider domain T> = 1536 x 1536, made possible thanks to the two-dimensional 
character of the model. A slight shift of the global threshold, from R\ ow ^170 for 256 x 128 
to « 171, was the price to be paid for the lowering of the spatial resolution by a factor of 
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Figure 1: Histograms of the values of E t (t) during 32 transients in the 32 x 32 system (a) 
and for the 32 sub-domains covering the 256 x 128 system (b). 

three (5x = 5z = 0.75, i.e., 12 x 6 collocation points per MFU) and of the temporal resolution 
by a factor of 20 (St = 0.2, i.e. five computations per turn-over time). This permitted data 
accumulation over long durations, typically of the order of 10 5 advective time units. 

Figure [2] (top) displays time series of the energy per unit surface E t contained in the 
perturbation and averaged over the whole domain T> for three experiments close to the 
global stability threshold. Experiments were performed by decreasing R from higher values 
to R = 172, then to R = 171 where a large laminar domain invaded the system, Fig. [2] lower 
left, and finally increasing R back to R = 171.5, for which all wide laminar domains receded 
leaving a whole distribution of small laminar patches, Fig. [2] lower right. The turbulent state 
then remained over a long period of time without any sign of decay. 

The statistical analysis of the flow patterns next proceeded with the detection of laminar 
domains. The local perturbation energy was coarse-grained over pixels of size 3x3 (i.e. 
1/2 MFU). A pixel was declared laminar when this coarse grained quantity was below 0.01 
and turbulent otherwise. The coarse graining procedure avoided the spurious detection of 
stagnation points in the velocity field. A Hoshen-Kopelman algorithm as implemented 
by Domany [2S] was used to identify the laminar clusters which were sorted according to 
their area measured in terms of pixels. Ranking clusters by decreasing size immediately gives 
cumulated distributions. The probability distributions shown in Figure [3] were obtained by 
differentiating the cumulated distributions. They were then slightly smoothed by resampling 
on bins of widths geometrically increasing with a factor 1.05. The pixel size, the "on/off" 
threshold and the smoothing factor resulted from a tedious trial-and-error procedure guided 
by the contradictory requirements of noise-level reduction and significance of the results. On 
the one hand, in contrast to the noise level, the slopes were not found sensitive to the value 
of the pixel-size/threshold combinations. On the other hand, the bin resizing factor was kept 
as close to one as possible in order to preserve the average value of the slopes visible in the 
raw data. 

Probability distribution of laminar cluster areas for R = 171.5 and R = 171 are displayed 
in Figure [3j The distribution for R = 172 resembles that for R = 171.5; extending a 
little less to large sizes, it is not presented for the sake of clarity. At higher values of R 
very few large laminar domains were detected and the distributions decayed faster, making 
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Figure 2: Top: time series of E t for the extra-wide domain T> = 1536 x 1536. Bottom: 
two b/w snap-shots of the solution for R = 171 (left) and R = 171.5 (right). Thresholding 
performed according to the rule given in the text. 
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Figure 3: Left: Probability distribution of cluster areas H(S) as a function of the surface S 
(numbers of laminar pixels). Right: S 3 U(S) as a function of S. R = 171.5: red/squares. 
R = 171: blue/circles. 



quantitative estimates unreliable. All these distributions were obtained by scanning a series 
of b/w snapshots such as those in Figure [2] bottom part. Since they were taken every At = 
3x 10 3 , the snapshots were sufficiently spaced in time to be statistically independent, whereas 
the large invading cluster correlating successive snapshots for R — 171 was systematically 
discarded from the statistics. All the distributions collapsed for laminar domains with areas 
below ten pixels (~ 5 MFUs). Tails with power law behavior and .R-dependent slopes were 
observed for the largest sizes. 

While the first momentum of the distribution of laminar domain sizes, (S) = 
SH(S)dS, remains defined, its second moment, (S 2 ) = S 2 U(S) dS, may diverge 
depending on whether or not the exponents a of the tails of the distributions are smaller or 
larger than 3. From the nucleation viewpoint, the divergence of the variance is in fact more 
meaningful than the divergence of the mean erroneously claimed in the version published 
initially [21] ■ Distributions of laminar domain sizes can be defined on time scales much longer 
than those of turbulent pattern renewal. If a laminar germ grows without bounds when it 
is larger than some critical size and recedes otherwise, then, as long as the variance remains 
somewhat smaller than the critical germ size, the probability of occurrence of a critical germ 
is negligible and turbulence remains sustained (R = 171.5, a ~ 3.4). When the variance 
diverge, a critical germ, whatever its size, will appear with finite probability and next will 
invade the system: turbulence is bound to decay for R = 171 (a ~ 2.8). The change of 
behavior is expected for 171.5 > R > 171 in our model and is conjectured to take place at 
R ~ 325 in Couette flow as reported by Bottin et al. [Hj, despite the reinterpretation of 
experimental data by Hof et al. [TT] . 

It might be argued that, despite its appealing properties, the model is not yet Navier- 
Stokes. In particular, due to its lack of cross-stream resolution, it cannot give informa- 
tion about the upper part of the transitional regime where oblique turbulent bands are 
observed [251126] . However it points toward a scenario that is likely to work in subcritical 
wall flows such as Ppf (quasi-lD), pCf (quasi-2D), or even the plane Poiseuille flow and the 
boundary layer flows (also quasi-2D) that all display regimes of laminar /turbulent coexis- 
tence at intermediate Reynolds numbers. STI is the process by which transient local chaos 



6 



is converted into sustained global spatiotemporal chaos. Two ingredients are necessary: a 
controllable decay rate of local states and a coupling with neighbors in physical space. Now, 
assume that, within the MFU assumption, the well-documented transient chaotic behavior 
has its lifetime primarily tuned by the Reynolds number. When they appear, nontrivial so- 
lutions have a short lifetime which leaves no time for contamination to work efficiently and 
so the flows relax to laminar. When R increases, the local deterministic dynamics in each 
MFU becomes more complicated and, as a result, the lifetime of transients increases. Due to 
the coupling, a genuine transition to sustained spatiotemporal chaos (weak turbulence) takes 
place at some given value R\ ow , whether or not the local lifetime, as obtained from the low- 
dimensional dynamical systems theory, diverges for R\ ow < R < oo. This transition would 
be within the class defined by STI and could be second-order (continuous) or first-order 
(discontinuous). Experiments tell us that, in the cases studied so far, it is first-order and 
thus deprived of any universality (correlation lengths remain finite at threshold). Abstract 
models considered up to now have been designed on lattices with local couplings and display 
both kinds of transitions in both one or two dimensions [20TT23] . The problem with this 
picture — solved by our (semi-)realistic Navier-Sokes modeling of pCf — is that couplings of 
hydrodynamic origin may have non-local effects linked to pressure, which makes size an im- 
portant issue. This should urge us to develop modeling further, not only in quasi-2D as was 
done here for pCf, but also in the quasi- ID case of Ppf by pushing further the simplification 
effort made in [30] in order to capture the spatiotemporal essence of the transition [32J. 
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